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Calculating Apparatus Having A Plurality of Stages 
Field of the invention: 

5 This invention relates to a clock distribution circuit for use with an apparatus 

having a large number of stages in what is often referred to as a pipeline arrangement. 

Background of the invention: 

10 It is becoming relatively common to exchange electronically stored documents 

between parties to a transaction, for instance via a widely distributed information network 
such as the Internet of the World Wide Web (WWW). A common problem with the 
^ Internet is a lack of secure communication channels. Thus, in order for hospitals, 
CO governments, banks, stockbrokers, and credit card companies to make use of the Internet, 

in 

|I 15 privacy and security must be ensured. One approach to solving the aforementioned 
It problem uses data encryption prior to transmission. In a prior art system, a host computer 
system is provided with an encryption unit, for example an encryption processor that is in 
electrical communication with at least a memory circuit for storing at least a private 
encryption key. When information is to be transmitted from the host computer system to 
%S 20 a recipient via the Internet and is of a confidential nature, the information is first passed 
to the encryption processor for encryption using the stored private key. Typically, a same 
private key is used every time a data encryption operation is performed. Alternatively, an 
encryption key is selected from a finite set of private encryption keys that is stored in the 
at least a memory circuit in electrical communication with the encryption processor. 

25 Of course, a data encryption operation that is performed by an encryption 

processor is a mathematical algorithm in which an input data value, for instance a hashed 
version of an electronic document, is the only variable value. It is, therefore, possible to 
optimize the encryption processor to perform a desired encryption function using a least 
amount of processor resources. Additionally, in the prior art encryption units the 
30 optimized encryption processor is typically separate from the microprocessor of the host 
computer system, because it is best optimized in this way. 
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Several standards exist today for privacy and strong authentication on the Internet 
through encryption/decryption. Typically, encryption/decryption is performed based on 
algorithms which are intended to allow data transfer over an open channel between 
5 parties while maintaining the privacy of the message contents. This is accomplished by 
encrypting the data using an encryption key by the sender and decrypting it using a 
decryption key by the receiver. In symmetric key cryptography, the encryption and 
decryption keys are the same. 

Encryption algorithms are typically classified into public-key and secret key 
10 algorithms. In secret-key algorithms, keys are secret whereas in public-key algorithms, 
one of the keys is known to the general public. Block ciphers are representative of the 
secret-key cryptosystems in use today. Usually, for block ciphers, symmetric keys are 

ij 

m used. A block cipher takes a block of data, typically 32-128 bits, as input data and 

m 

produces the same number of bits as output data. The encryption and decryption 
E ? ™15 operations are performed using the key, having a length typically in the range of 56-128 
£n bits. The encryption algorithm is designed such that it is very difficult to decrypt a 
~ message without knowing the key. 

In In addition to block ciphers, Internet security protocols also rely on public-key 

*fj based algorithms. A public key cryptosystem such as the Rivest, Shamir, Adelman 
p20 (RSA) cryptosystem described in U.S. Pat. No. 5,144,667 issued to Pogue and Rivest 
uses two keys, one of which is secret - private - and the other of which is publicly 
available. Once someone publishes a public key, anyone may send that person a secret 
message encrypted using that public key; however, decryption of the message can only be 
accomplished by use of the private key. The advantage of such public-key encryption is 
25 private keys are not distributed to all parties of a conversation beforehand. In contrast, 
when symmetric encryption is used, multiple secret keys are generated, one for each 
party intended to receive a message, and each secret key is privately communicated. 
Attempting to distribute secret keys in a secure fashion results in a similar problem as 
that faced in sending the message using only secret-key encryption; this is typically 
30 referred to as the key distribution problem. 
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Key exchange is another application of public-key techniques. In a key exchange 
protocol, two parties can agree on a secret key even if their conversation is intercepted by 
a third party. The Diffie-Hellman exponential key exchange method, described in U.S. 
Pat. No. 4,200,770, is an example of such a protocol. 

5 Most public-key algorithms, such as RSA and Diffie-Hellman key exchange, are 

based on modular exponentiation, which is the computation of a x mod p. This 
expression means "multiply a by itself x times, divide the answer by p, and take the 
remainder." This is very computationally expensive to perform, for the following reason. 
In order to perform this operation, many repeated multiplication operations and division 
10 operations are required. Techniques such as Montgomery's method, described in 

"Modular Multiplication Without Trial Division," from Mathematics of Computation, 
p Vol. 44, No. 170 of April 1985, can reduce the number of division operations required 

^ but do not overcome this overall computational expense. In addition, for present day 

CO 

in encryption systems the numbers used are very large (typically 1024 bits or more), so the 

il 15 multiply and divide instructions found in common CPUs cannot be used directly. 

^ Instead, special algorithms that break down the large multiplication operations and 

-J 

s division operations into operations small enough to be performed on a CPU are used. 

ti These algorithms usually have a run time proportional to the square of the number of 

O machine words involved. These factors result in multiplication of large numbers being a 

S3. 

□ 20 very slow operation. For example, a Pentium® processor can perform a 32x32-bit 

multiply in 10 clock cycles. A 2048-bit number can be represented in 64 32-bit words. A 
2048x2048-bit multiply requires 64x64 separate 32x32-bit multiplication operations, 
which takes 40960 clocks on the Pentium® processor. An exponentiation with a 2048-bit 
exponent requires up to 4096 multiplication operations if done in the straightforward 
25 fashion, which requires about 167 million clock cycles. If the Pentium processor is 
running at 166 MHZ, the entire operation requires roughly one second. Of course, the 
division operations add further time to the overall computation times. Clearly, a common 
CPU such as a Pentium cannot expect to do key generation and exchange at any great 
rate. 
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Pipeline processors comprising a plurality of separate processing elements 
arranged in a serial array, and in particular a large number of processing elements, are 
known in the prior art and are particularly well suited for executing data encryption 
algorithms. Two types of pipeline processor are known: processors of an in-one-end- 
5 and-out-the-other nature, wherein there is a single processing direction; and, bi- 
directional processors of an in-and-out-the-same-end nature, wherein there is a forward 
processing direction and a return processing direction. Considering a specific example of 
a bi-directional pipeline processor, a first data block is read from a memory buffer into a 
first processing element of the serial array, which element performs a first stage of 
10 processing and then passes the first data block on to a second processing element. The 
second processing element performs a second stage of processing while, in parallel, the 
first processing element reads a second data block from the memory buffer and performs 

□ a same first processing stage on the second data block. In turn, each data block 

propagates in a step-by-step fashion from one processing element to a next processing 

Fl5 element along the forward processing direction of the serial array. At each step, there is a 
processing stage that performs a same mathematical operation on each data block that is 

JpJ provided thereto. Simultaneously, a result that is calculated at each processing element is 
provided to a previous processing element of the serial array, with respect to the return 

in processing direction, which results comprise in aggregate the processed data returned by 

*j&0 the encryption processor. This assembly-line approach to data processing, using a large 

ted" 

Q number of processing elements, is a very efficient way of performing the computationally 
expensive data encryption algorithms described previously. Of course, the application of 
pipeline processors for performing computationally expensive processing operations is 
other than limited strictly to data encryption algorithms, which have been discussed in 
25 detail only by way of example. 



It is a disadvantage of the prior art bi-directional pipeline processors that each 
processing element of a serial array must be time-synchronized with every other 
processing element of a same serial array. Time-synchronization between processing 
30 elements is necessary for the control of timing the gating of data blocks from one 

processor element to a next processor element in the forward direction, and for timing the 
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gating of processed data from one processor element to a previous processor element in 
the return direction. A clock typically controls the progression of data blocks along the 
pipeline in each one of the forward direction and the return direction. Unfortunately 
without careful clock distribution design, as a clock signal progresses along the pipeline 
there are incremental delays between each stage, as for example delays caused by the 
resistance and capacitance that is inherent in the clock circuit. In earlier, slower acting 
pipeline processors, such delays were not important, and did not adversely affect the 
overall operation, or calculation. With faster operation, these delays are becoming 
significant, requiring more accurate and precise clock distribution methods. 



Further, in order to read data from a memory buffer, for example data for 
processing by the pipeline processor, the first processing stage in the serial array must 
? 3 also be time-synchronized with the memory buffer. This further encourages 
synchronous clock distribution within a pipeline processor. 



It would be advantageous to provide a system and a method for processing data 
using a pipeline processor absent a need to synchronize a distributed clock value that is 
provided to each processing element of the pipeline processor. Such a system would be 
easily implemented using a relatively simple circuit design, in which large blocks of 



, R 20 processor elements are fabricated from a series of processor element sub-units. 



Object of the Invention 

It is an object of the invention to provide a pipeline processor absent a 
synchronous clock signal for all processing elements. 

Summary of the invention 



In its broadest concept, the invention provides a calculating apparatus having a 
plurality of stages in an extended pipeline array, arranged in a series of side-by-side sub- 
arrays, and a clock conductor extending in a sinuous form alongside the array, connected 
30 to each stage. The array can be in the form of sections each having input and output 



5 



Doc. No. 4747 US 




Patent 



access whereby the whole array or sections of the array can process data. The apparatus 
has forward and return paths and can be arranged so that the shortest calculation taking 
place in a stage is arranged to take place in the return path. 

5 In accordance with another embodiment of the invention there is provided an 

apparatus for processing data comprising: 

a plurality of individual processing elements arranged in a serial array wherein a 
first processing element precedes a second processing element which precedes an nth 
processing element; and, 
10 a clock distribution circuit in electrical communication with each processing 

element of the plurality of individual processing elements in the serial array such that, in 
use, a clock signal propagated along the clock distribution circuit arrives at each 
E» processing element delayed relative to the clock signal arriving at a preceding processing 

£Q element; 

In 

j7 15 wherein a time equal to an exact number of clock cycles, k, where k is 

H greater than zero, from when the data is clocked into a processing element to when the 
•Jj data is clocked in by a subsequent processing element is insufficient for providing 
L accurate output data from the processing element but wherein the same time with the 
111 additional delay is sufficient and wherein new data to be processed is clocked in by the 
same processing element after the exact number of clock cycles, k. 

C3 

In accordance with another embodiment of the invention there is provided a 
switchable processing element comprising: 

a first port for receiving a first clock signal; 
25 a second port for receiving a second other clock signal; 

a switch operable between two modes for selecting one of the first clock signal 
and the second other clock signal; and 

wherein the selected one of the first clock signal and the second other clock signal 
is provided to the processing element. 

30 
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In accordance with another aspect of the invention there is provided a method for 
processing data comprising the steps of: 

(a) providing a pipeline processor including a plurality of individual 
processing elements arranged in a serial array such that a first processing element 

5 precedes a second processing element which precedes an nth processing element; 

(b) providing a clock signal to each processing element of the plurality of 
individual processing elements in the serial array such that the clock signal arrives at each 
individual processing element beyond the first processing element delayed relative to the 
clock signal arriving at a preceding processing element; 

10 (c) providing data to the first processing element for processing therein; and, 

(d) propagating the data to at least a next processing element for additional 
processing therein, 

CD wherein the clock signal provided to an element in the plurality of 

HJ individual processing elements is delayed relative to the clock signal provided to another 

if! 

r":15 element of the plurality of individual processing elements by a substantial amount 
H relative to the clock period. 

2 In accordance with another embodiment of the invention there is provided a 

0 

If! method for processing data within a pipeline processor comprising the steps of: 

n 

*i20 (a) providing a clock signal in a first direction along a first portion of the 

O pipeline processor having a number, n, processing elements such that the clock signal 
arrives at each individual processing element beyond the first processing element of the 
first portion delayed relative to the clock signal arriving at a preceding processing 
element of the same first portion; 
25 (b) providing a clock signal in a second substantially opposite direction along 

a second other portion of the pipeline processor having a same number, n, processing 
elements such that the clock signal arrives at each individual processing element beyond 
the first processing element of the second other portion delayed relative to the clock 
signal arriving at a preceding processing element of the same second other portion; 
30 (c) providing data to the first processing element of the first portion of the 

pipeline processor for processing therein; 
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wherein the delay to the last processing element of the first portion is an 
approximately same delay as the delay to the last processing element of the second 
portion, such that at center of the pipeline processor the two adjacent processing elements 
are in synchronization. 

5 

In accordance with yet another aspect of the invention there is provided a macro 
for use in layout of an apparatus for processing data comprising: 

a plurality of individual processing elements arranged serially and having a clock 
10 input conductor and a clock output conductor, the clock input conductor in 

communication with a clock conductor having increased length from the clock input 
conductor to each subsequent element within the within the plurality of individual 
□ processing elements and wherein the clock conductor has decreased length from the clock 

'it output conductor to each subsequent element within the within the plurality of individual 

ly 

Ifl 15 processing elements, 

^ wherein the clock input conductor and output conductor are arranged such that 

adjacently placed macros form space efficient blocks within a layout and such that the 

s input clock conductor of one macro and the out clock conductor of an adjacent macro 

?S when coupled have approximately a same conductor path length as the conductor path 

*3 20 length between adjacent elements within a same macro when the macros are disposed in a 

Sj.ijj 

p predetermined space efficient placement. 

Brief description of the drawings 

25 The invention will be readily understood by the following description of preferred 

embodiments, in conjunction with the following drawings, in which: 
Figure 1 shows a simplified block diagram of a first preferred embodiment of a pipeline 
processor according to the present invention; 

Figure 2 shows a simplified block diagram of an array of processor elements in electrical 
30 communication with a distributed clock circuit according to the present invention; 
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Figure 3 shows a timing diagram for gating information to a plurality of processor 
elements in a prior art pipeline processor; 

Figure 4 shows a timing diagram for gating information to a plurality of processor 
elements in a pipeline processor, according to the present invention; 
5 Figure 5 shows individual timing diagrams for three adjacent processor elements within a 
same processor array according to the present invention; 

Figure 6 shows a simplified block diagram of a second preferred embodiment of a 
pipeline processor according to the present invention; 

Figure 7 shows a simplified block diagram of a third preferred embodiment of a pipeline 
10 processor according to the present invention; 

Figure 8a shows a simplified block diagram of a processor element having a clock 
switching circuit and operating in a first mode according to the present invention; 
Figure 8b shows a simplified block diagram of a processor element having a clock 
y switching circuit and operating in a second mode according to the present invention 

in 

M> 15 Fig. 9 is a simplified block diagram of macro blocks of processor units arranged for 

I* 

ff. providing a snaking clock signal from unit to unit; 

Figure 10 is a block diagram of a resource efficient processing element design for use in a 

C3 pipeline array processor for performing encryption functions; 

in 

n Figure 1 1 is a block diagram of a systolic array for modular multiplication; 
20 Figure 12 is a block diagram of a single unit with its input pathways shown; 

□ 

^ Figure 13 is a block diagram of a DP RAM Z unit; 
Figure 14 is a block diagram of an Exp RAM unit; 
Figure 15 is a block diagram of a Prec RAM unit; 

Figure 16 is a block diagram of a speed efficient processing element design for use in a 
25 pipeline array processor for performing encryption functions; 

Figure 17 is a block diagram of a systolic array for modular multiplication; 
Figure 18 is a block diagram of a single unit with its input pathways shown; and, 
Figure 19 is a block diagram of a DP RAM Z unit. 
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Detailed description of the invention 



10 



The present invention is concerned with the reduction of time delays between 
stages. The result is obtained by positioning a clock conductor in the proximity of the 
various stages, as by snaking the conductor alongside the stages. Thus the clock delay is 
now substantially small between adjacent elements without a need for proper inter- 
element synchronization. A further advantage is realized when a consistent time delay is 
provided between adjacent elements in that interconnection between stages other than 
those immediately adjacent is possible. 

A further advantage is that, if desired, instead of the entire array of stages being 
used for a large calculation, the array can be subdivided, for example into halves or 

JpL quarters, such that more than one calculation is carried out at a same time. 

ff\ 

in 

U 15 Referring to Figure 1, shown is a simplified block diagram of a pipeline processor 

\t 7 in electrical communication with a real time clock 1 via a hardware connection 2, 

In 

*P according to a first embodiment of the present invention. The pipeline processor 7 
h\ includes a plurality of arrays 4a, 4b and 5 of processor elements (processor elements not 
shown), for instance, arrays 4a and 4b each has 256 processing elements and array 5 has 
vO 20 512 processing elements. An input/output port 9 is separately in communication with the 
first processing element of each array 4a, 4b and 5, for receiving data for processing by 
the pipeline processor 7, for example from a client station (not shown) that is also in 
operative communication with the port 9. A clock conductor 3, in electrical 
communication with clock source 1 via hardware connection 2, is provided in the form of 
25 a distributed clock circuit extending in a sinuous form alongside each of arrays 4a, 4b 
and 5. The clock conductor 3 is also separately in electrical communication with each 
individual processor element of the arrays 4a, 4b and 5. 

Referring to Figure 2, shown is a simplified block diagram of a serial array of 
30 processor elements 8 1 , 8 2 , 8 3 , ... , 8"* 1 and 8", the individual processor elements 8 

comprising in aggregate the array 4a of pipeline processor 7 in Figure 1. Each processor 
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element 8 is separately in electrical communication with the clock conductor 3 via a 
connection 10. The clock conductor 3 is also in electrical communication with a clock 
generator circuit, the clock source, via hardware connection 2. An input/output port 9 in 
communication with the first processing element of array 4a is for receiving data 
provided by a client station (not shown), also in operative communication with 
input/output port 9, the data for processing by the array 4a. 

In operation, data is provided by the client station at port 9, for example as a 
stream of individual blocks of data which comprise in aggregate a complete data file. 
The first processor element 8 ! in array 4a receives a first data block via port 9 and 
performs a predetermined first processing stage thereon. Of course, first processor 
element 8 l is time-synchronized with a memory buffer (not shown) of port 9 such that the 
stream of data blocks is gated to first processor element 8 1 in synchronization. For 
example, clock conductor 3 provides a time signal from real time clock 1, the time signal 
arriving at first processor element 8 1 at a predetermined time relative to a clock signal of 
the memory buffer. At the end of a first processing cycle, first processor element 8 1 
receives a second data block via port 9. At a same time the first processing element 8 1 
provides an output from the first data block along a forward processing-path to second 
processor element 8 2 . Additionally, the first processor element 8 1 provides a second result 
calculated therein along a return processing-path to the buffer of port 9. 

During a second processing cycle, first processor element 8 1 performs a same first 
processing operation on the second data block and second processor element 8 performs 
a second processing operation on the first data block. At the end of the second 
processing cycle, the result of processing on the first data block is propagated along the 

2 3 

forward processing path between the second and the third processor elements 8 and 8 , 
respectively. Simultaneously, the results of processing of the second data block is 
propagated along the forward processing path between the first and the second processor 
elements 8 1 and 8 2 , respectively. Additionally, the second processor element 8 2 provides 
a result calculated therein along a return processing-path to the first processor element 8 1 . 
Of course, simultaneously gating data blocks along the forward processing-path and 



11 



Doc. No. 47-17 US 




Patent 



along the return processing-path between adjacent processor elements requires 
synchronous timing. For instance, it is critical that the processing operations that are 
performed along both processing-paths are complete prior to the data being propagated in 
either direction. 

Referring to Figure 3, shown is timing diagram for gating information to a 
plurality of processor elements in a prior art pipeline processor. By way of example, 
individual timing diagrams for a first five processor elements, denoted 1, 2, 3, 4 and 5, 
respectively, are shown. Each clock cycle is denoted by a pair of letters, for example AB, 
CD, EF, etc. It is assumed for the purpose of this description that information is gated to 
and from each processor element at a "rising edge" of any clock cycle. For instance, 
along the forward processing path processor element 1 gates in a first block of data at 
"rising edge" AB and processes the first block of data during one complete clock cycle. 
Similarly, processor element 2 gates in the first block of data from processing element 1 
at "rising edge" CD and processes the first block of data during one complete clock cycle. 
Additionally, along the return processing-path, processor element 1 gates in a block of 
processed data prom processor element 2 at "rising edge" EF. 

Of course, the clock cycle rate of the prior art system is at least as long as the 
longest processing time required at each stage along one of the forward and the return 
processing paths. For example, a data stream propagates along the serial array in a 
stepwise fashion, and processing must be completed at every step before the data can be 
propagated again. Thus if processing occurs in a shorter period of time along the return 
processing path compared to the forward processing path, then a delay is introduced at 
every stage along the reverse processing path in order to allow the processing to be 
completed along the forward processing path. 

Additionally, as is apparent from Figure 3, every processor element must be 
synchronized with every other processor element of the array. For instance the clock 1 of 
Figure 1 must be distributed everywhere along the array in phase. This typically is a 
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complex problem that is costly and difficult to solve. The solutions are usually a hybrid 
of hardware design and integrated circuit topology design and analysis. 

An approach to overcoming the problem of clock distribution is a technique 
5 wherein a first processor provides a clock signal to a second processor and from there it is 
provided to a third processor and so forth. Thus, between adjacent elements, 
synchronization exists but, between distant elements, synchronization is not assured. 
Unfortunately, this method of avoiding clock synchronization is performed absent a 
global clock and, as such, a clock is passed between every two elements requiring data 
10 communication therebetween resulting in a different clock distribution problem. 

Referring to Figure 4, shown is a timing diagram for gating information to a 

^3 plurality of processor elements in a pipeline processor, according to the present invention. 

CO By way of example, the individual timing diagrams for a subset of a serial array 

jj\ 15 comprising the first ten processor elements, denoted 1, 2, 3, 4, 5, 6, 7, 8, 9, and 10, 

^ respectively, are shown. Each clock cycle is denoted by a pair of letters, for example AB, 

kg CD, EF, etc. It is assumed for the purpose of this discussion that information is gated 

^ into and out of each processor element at a "rising edge" of a clock cycle. For instance, 

.lor 

in along the forward processing path processor element 1 gates in a first block of data at 

O 

In 20 "rising edge" AB and processes the first block of data during one complete clock cycle. 
W Similarly, processor element 2 gates in the first block of data from processing element 1 

teas 

at "rising edge" CD and processes the first block of data during one complete clock cycle. 
Additionally, along the return processing-path, processor element 1 gates in a block of 
processed data prom processor element 2 at "rising edge" EF. It is further assumed for 
25 the purpose of this discussion that the processing operation requiring the greatest amount 
of time to be completed at any processor element is along the forward processing-path. 
Of course, as indicated by the diagonal lines in Fig. 4, the rising edge AB occurs at 
different times for different processing elements. 

30 Referring still to Figure 4, each timing diagram is offset slightly from the timing 

diagram for a previous processor element by an amount, 5, equal to an incremental delay 
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of the clock signal reaching that processing element. Due to capacitance and resistance 
that is inherent in the circuitry comprising the clock conductor, the finite period of time, 
5, elapses between the arrival of the time signal at the first processor element and the 
arrival of the time signal at the second processor element. Alternatively, the clock is 
5 intentionally delayed between provision to different processing elements. Thus, the time- 
synchronization between processor element 1 and processor element 2 is offset by the 
amount 6. Similarly, the time-synchronization between each of the remaining pairs of 
adjacent processor elements also is offset, for example by a same amount 8. 
Alternatively, the offset amount is different but within known tolerances. 

10 

Still referring to Figure 4, the individual clock cycles are shorter than the clock 
C3 cycles of the prior art timing diagrams shown in Figure 3 for a same processing 

(Q operation. This would seem to imply that there is insufficient time for the processor 

i ^ 

f ; elements to complete the processing operations along the forward processing-path prior 

f£ 

15 to gating in new data. For example, in Figure 3 the clock cycle is at least as long as the 
Jpj longest processing operation, which operation is arranged to occur along the forward 
5 =i path. In the present embodiment, however, there is an incrementally increasing delay of 
|fl the arrival of the clock signal at each processing element beyond processor element 1. In 
^ effect, this delay provides additional time for processing to be completed at, for example, 
0 20 processor element 2 in a forward processing path before a next block of data is gated in at 
processing block 3 from processor element 2. Advantageously, the minimum length of 
an individual clock cycle is reduced to a length of time equal to the time required to 
complete the longest processing operation less the length of the clock delay between 
elements in the path requiring longer processing times - here the forward path. Then, 
25 along the forward processing path more than one full clock cycle elapses between gating 
a block of data into a processor element and gating the processed block of data from that 
processor element into a next processor element. Further, along the return processing 
path less than one full clock cycle elapses between gating a block of data into a processor 
element and gating the processed block of data into a next processor element (previous in 
30 the forward path). The invention provides what can be termed "catch up" in the return 
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processing-path. Thus, the overall cycle time is less than the time required in one 
direction of processing but at least an average of the processing time required in each of 
the two directions. 

5 Referring to Figure 5, shown are three individual timing diagrams for three 

adjacent processor elements, denoted 3, 4 and 5, according to the present invention. A 
first data block is gated into processor element 4 at 100 and is processed by processor 
element 4 during clock cycle FG. For example, processor element 4 reads the first data 
block from an output port of processor element 3, the first data block having been gated 
10 into processor element 3 at 101. Processor element 4 also makes the first data block 

available to processor element 5, for example processor element 4 provides the first data 
block to an output port thereof and the first data block is read by processor element 5 at 
W 104. Clearly, steps 101, 100 and 104 comprise a portion of the forward processing-path. 
Cy As is obvious from Figure 5, a period of time that is longer than one complete clock cycle 
^ 15 elapses between gating a block of data into a processor element and gating a block of data 

S' ■ 

^ resulting from processing of the same block of data into a next processor element along 
CP 

v3 the forward processing-path. 

a 

C3 

in Similarly, the steps 102, 100 and 103 comprise a portion of the reverse 

20 processing-path, wherein a data block including data processed by a processor element is 

^ provided to a previous processor element of the array. As is obvious from Figure 5, a 
period of time that is shorter than one complete clock cycle elapses between gating a 
processed block of data into a processor element and gating the further processed block 
of data into a next processor element along the return processing-path. Advantageously, 
25 the processing delay that accumulates along the forward processing-path is "caught-up" 
along the return processing-path. This is a phenomenon that is referred to as "bi- 
directional averaging". Further, since the length of the clock cycle time is reduced in the 
present invention, an overall advantage in increased processing speed over prior art bi- 
directional pipeline processors is realized. 

30 
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It is an advantage of the present invention that each processor element needs only 
to communicate with two adjacent elements, such that an exact delay is always 
determinable and can easily be maintained within predetermined limits. It is a further 
advantage of the present invention that it is possible to isolate the circuit design to n 
5 adjacent processor elements, such that the entire pipeline processor is fabricated by 
laying down a series of n element "macros". Of course, every once in a while it is 
necessary to connect one macro block to another, requiring additional circuitry to cope 
with an extra delay between processor elements of different macro blocks. Alternatively, 
macros are designed for ease of interconnection such that a macro begins and ends in a 
10 fashion compatible with positioning another identical macro adjacent thereto for 
continued similar performance. In Fig. 9, a diagram of 2 macro blocks 91 and 92 
according to the invention is shown. The macro blocks can be arranged in any of a series 
S3 of arrangements as shown providing approximately consistent pathway delays between 

f Q processing elements. 

m 

7715 

M Referring to Figure 6, shown is a simplified block diagram of a pipeline processor 

Co 

^9 12 according to a second preferred embodiment of the present invention. The pipeline 
L processor 12 includes a plurality of arrays 4a, 4b and 5 of processor elements (processor 

%U elements not shown), for instance, arrays 4a and 4b each having 256 processing elements 

t i 

Ifi 20 and array 5 having 512 processing elements. Dotted lines 6a and 6b indicate optional 
^ electrical coupling for providing electrical communication between the 256 th processing 
element of array 4a and the 256 th element of array 4b, and between the 1 st element of 
array 4b and the 1 st element of array 5, respectively. A distributed clock circuit 3 is 
separately in electrical communication with each processor element of the arrays 4a, 4b 
25 and 5. Also shown in Figure 6 is a clock generator 1 in electrical communication with 
pipeline processor 12 via a hardware connection 2. An input/output port 9 in 
communication with the first processing element of each array 4a, 4b, and 5 is for 
receiving data provided by a client station (not shown), also in operative communication 
with input/output port 9, the data for processing by an indicated one of the arrays 4a, 4b, 
30 and 5. 
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Referring to Figure 7, shown is a simplified block diagram of a pipeline processor 
13 according to a third preferred embodiment of the present invention. The pipeline 
processor 13 includes a plurality of arrays 4a, 4b and 5 of processor elements (processor 
elements not shown), for instance, arrays 4a and 4b each having 256 processing elements 
5 and array 5 having 512 processing elements. The 256 th processing element of array 4a 
and the 256 th element of array 4b are in electrical communication via the hardware 
connection 11a, and the 1 st element of array 4b and the 1 st element of array 5 are in 
electrical communication via the hardware connection lib, respectively. A distributed 
clock circuit 3 is separately in electrical communication with each processor element (not 
10 shown) of the arrays 4a, 4b and 5. Also shown in Figure 7 is a real time clock 1 in 
electrical communication with pipeline processor 13 via a hardware connection 2. An 
input/output port 9 in communication with the first processing element of array 4a is for 
G receiving data provided by a client station (not shown), also in operative communication 

m with input/output port 9, the data for processing by the serial arrangement of the arrays 

in 

15 4a, 4b, and 5. Optionally, separate inputs (not shown) are provided for gating data 

H directly to at least a processor element other than the 1 st element of array 4a. 

m 
^0 



o 



The pipeline processors 12 and 13 of Figures 6 and 7, respectively, are operable in 

10 mode wherein data gated into the 256 th processor element of the array 4a is made 

O 

In 20 available to the 256 processor element of array 4b. For instance, when more than 256 
l : J processor elements are required for a particular processing operation, the effective length 
of the processor array is increased by continuing the processing operation within a second 
different array. Of course, when more than 512 processor elements are required for a 
particular processing operation, the effective length of the processor array is increased by 
25 continuing the processing operation within a third different array. For example, either 
one of the pipeline processors shown in Figures 6 and 7 are operable for performing: 256 
bit encryption using a single array; 512 bit encryption using two different arrays; and, 
1024 bit encryption using all three different arrays. Of course, optionally the 256 th 
processor element of array 4a is coupled to the 1 st element of array 4b, but then both the 
30 256 th element of array 4a and the 1 st element of array 4b must be synchronized with each 
other and with the buffer. Such synchronization requirements increase the circuit design 
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complexity due to the critical need for a uniform distributed clock. Also, in most pipeline 
processor arrangements it is necessary that each element provide processing operations 
during each clock cycle and often, clock synchronization imposes a wait state which 
would cause the 257 th element in the array to process data one clock cycle later than the 
earlier elements. 

Of course, when the 256 th element of array 4a is coupled to the 256 th element of 
array 4b, either optionally as shown in Figure 6 or permanently as shown in Figure 7, the 
advantage of "bi-directional averaging" is lost. Advantageously, however, a plurality of 
separate arrays of processor elements, each array preferably comprising a same number 
of processor elements, is connectable in such a head-to-tail fashion. Then, the clock 
signal is delayed progressively along every second array, but catches-up again in 
between. 

Of course, since clock distribution is not a significant concern and delays in clock 
distribution are well supported, the clock signal is optionally switched into each 
processing element such that the clock is provided from one of two clocking sources. 
Then, with a processor circuit configuration similar to that of Fig. 7, the clock is switched 
in direction for the second processor array and provided through coupling 11a. Thus the 
advantages of "catch up" are maintained and synchronization between adjacent arrays is 
obviated. Further, such a configuration supports arrays of various length that are 
couplable one to another to form longer arrays when needed without a necessity for clock 
synchronization therebetween. Here, every processing element within the second array 
requires two clock sources - one from a preceding element in a first direction and another 
from a preceding element in a second other direction. Since clocks are delayed between 
processing elements, the switching circuit merely acts to impart a portion or all of the 
necessary delay to the clock signal. 

Referring to Fig. 8, a processing element is shown having a clock switching 
circuit for use according to the present embodiment. A first clock signal is provided at 
port 81. A second other clock signal is provided at port 82. Since, in use, the clock only 
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propagates along one direction, the ports 81 and 82 are optionally bi-directional ports. 
Each port is coupled to a clock driver 84 and 83 respectively. The ports are also coupled 
to a switch 85 for providing only one selected clock along a clock conductor 86 to the 
processing element 87. The clock is also provided to the two drivers only one of which is 
5 enabled. In this way, each element works to propagate a clock signal in one direction 
selectable from two available directions of clock propagation. 

Advantageously, since it is known when a processor will complete processing, it 
becomes possible to allocate that processor to processing downstream of another 
10 processor. For example, assuming the processor 4a has processing elements for 

processing 256 bit operations and begins processing a 256 bit operation. Assume 4b is a 
similar processor. If, sometime after processing element 4a commences processing and 
before it is completed a processing request for a 512 bit operation arrives, it is possible to 
begin the operation on processing array 4b knowing that by the time data has propagated 
15 to the last element of processing array 4a, that element will have completed processing of 
the processing job in current processing. This improves overall system performance by 
reducing downtime of a processor while awaiting other processors to be available to 
JL support concatenated array processing. 

in 

j/g 20 Montgomery based Pipeline processing of encryption data 



£3 
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Applying Montgomery's algorithm, the cost of a modular exponentiation is 
reduced to a series of additions of very long integers. To avoid carry propagation in 
multiplication/addition architectures several solutions are known. These use 
25 Montgomery's algorithm, in combination with a redundant radix number system or a 
Residue Number System. 

In S.E.Eldridge and C.D.Walter.Hardware implementation of Montgomery's 
modular multiplication algorithm. IEEE Transactions on Computers, 42(6):693-699,July 
1993, Montgomery's modular multiplication algorithm is adapted for an efficient 
30 hardware implementation. A gain in speed results from a higher clock frequency, due to 
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simpler combinatorial logic. Compared to previous techniques based on Brickell's 
Algorithm, a speed-up factor of two was reported. 

The Research Laboratory of Digital Equipment Corp. reported in J. E. Vuillemin, 
P. Bertin, D. Roncin, M. Shand, H.H. Touati, and P. Boucard. Programmable active 
memories: Reconfigurable systems come of age. IEEE Transactions on VLSI Systems, 
4(1): 56-69, March 1996 and M. Shand and J. Vuillemin. Fast implementations of RSA 
cryptography. In Proceedings 11th IEEE Symposium on Computer Arithmetic, pages 
252-259, 1993, an array of 16 XILINX 3090 FPGAs using several speed-up methods 
including the Chinese remainder theorem, asynchronous carry completion adder, and a 
windowing exponentiation method is used to implement modular exponentiation. The 
implementation computes a 970bit RSA decryption at a rate of 185kb/s (5.2ms per 970 
bit decryption) and a 512 bit RSA decryption in excess of 300 kb/s (1.7ms per 512 bit 
decryption). A drawback of this solution is that the binary representation of the modulus 
is hardwired into the logic representation so that the architecture must be reconfigured 
with every new modulus. 

The problem of using high radices in Montgomery's modular multiplication 
algorithm is a more complex determination of a quotient. This behavior renders a 
pipelined execution of the algorithm other than straightforward. In H.Orup. Simplifying 
quotient determination in high-radix modular multiplication. In Proceedings 12th 
Symposium on Computer Arithmetic, pages 193-9, 1995, the algorithm is rewritten to 
avoid any operation involved in the quotient determination. The necessary pre- 
computation is performed only once for a given modulus. 

P. A. Wang in the article New VLSI architectures of RSA public key crypto 
systems. In Proceedings of 1997 IEEE International Symposium on Circuits and Systems, 
volume 3, pages 2040-3, 1997 proposes a novel VLSI architecture for Montgomery's 
modular multiplication algorithm. The critical path that determines the clock speed is 
pipelined. This is done by interleaving each iteration of the algorithm. Compared to 
previous propositions, an improvement of the time-area product of a factor two was 
reported. 
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J.Bajard, L.Didier, and P.Kornerup in the article An RNS Montgomery modular 
multiplication algorithm. IEEE Transactions on Computers, 47(7) : 766 -76, July 1998, 
describe a new approach using a Residue Number System (RNS). The algorithm is 
implemented with n moduli in the RNS on n reasonably simple processors. The resulting 
processing time is O(n). 

Of course, most of the references cited above relate to hardware implementations 
of processors that have little or no flexibility. 

There have also been a number of proposals for systolic array architectures for 
modular arithmetic. These vary in terms of complexity and flexibility. 

In E. F. Brickell. A survey of hardware implementations of RSA. In Advances in 
Cryptology —CRYPTO '89, pages 368-70.Springer-Verlag, 1990, E.F. Brickell 
summarizes the chips available in 1990 for performing RSA encryption. In 

In N. Takagi. A radix-4 modular multiplication hardware algorithm efficient for 
iterative modular multiplication operations . In Proceedings 10th IEEE Symposium on 
Computer Arithmetic, pages 35-42, 1991, the author proposes a radix-^ hardware 
algorithm. A redundant number representation is used and the propagation of carries in 
additions is therefore avoided. A processing speed-up of about six times compared to 
previous work is reported. 

More recently an approach has been presented that utilizes pre-computed 
complements of the modulus and is based on the iterative Horner's rule in J. Yong-Yin 
and W. P. Burleson. VLSI array algorithms and architectures for RSA modular 
multiplication. IEEE Transactions on VLSI Systems, 5(2): 211-17, Jun 1997. Compared 
to Montgomery's algorithms these approaches use the most significant bits of an 
intermediate result to decide which multiples of the modulus to subtract. The drawback of 
these solutions is that they either need a large amount of storage space or many clock 
cycles to complete a modular multiplication. 

The most popular algorithm for modular exponentiation is the square & multiply 
algorithm. Public-key encryption systems are, typically, based on modular 
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exponentiation or repeated point addition. Both operations are in their most basic forms 
done by the square and multiply algorithm. 

Method 1.1 compute Z = X E mod M, where E = X"^2' f e r . e {0,l} 

1. Z = X 

2. FOR i = n - 2 down to 0 DO 

3. Z = Z 2 mod M 

4. IF ei = 1 THEN Z = Z • X mod M 

5. END FOR 

Method 1.1 takes 2(n-l) operations in the worst case and 1.5(n-l) on average. To 
compute a squaring and a multiplication in parallel, the following version of the square & 
multiply method can be used: 

Method 1.2 computes P = X E mod M, where E = 2' , e. e {0,l} 

1. P 0 = l f Zo = X 

2. FOR i = 0 to n - 1 DO 

3. Z I+1 =Zi 2 modM 

4. IF ^ = 1 THEN P i+1 = Pi • Z x mod M 
ELSE P i+ i = Pi 

5. END FOR 

Method 1.2 takes 2n operations in the worst case and 1.5n on average. A speed- 
up is achieved by applying the 1 - ary method, such as that disclosed in D. E. Knuth, The 
Art of Computer Programming. Volume 2: Seminumerical Algorithms. Addison- Wesley, 
Reading, Massachusetts, 2nd edition, 1981, which is a generalization of Method 1.1. The 
1 - ary method processes 1 exponent bits at a time. The drawback here is that (2 1 - 2) 
multiples of X must be pre-computed and stored. A reduction to 2 1 " 1 pre-computations is 
possible. The resulting complexity is roughly n/1 multiplication operations and n 
squaring operations. 
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As shown above, modular exponentiation is reduced to a series of modular 
multiplication operations and squaring steps using the Montgomery method. The method 
for modular multiplication described below was proposed by P. L. Montgomery in P. L. 
Montgomery. Modular multiplication without trial division. Mathematics of 
5 Computation, 44(170): 519-21, April 1985. It is a method for multiplying two integers 
modulo M, while avoiding division by M. The idea is to transform the integers in m- 
residues and compute the multiplication with these m-residues. In the end, the 
representations are transformed back to a normal representation thereof. This approach is 
only beneficial when a series of multiplication operations in the transform domain are 
10 computed (e.g., modular exponentiation). 

To compute the Montgomery multiplication, a radix R > M, with gcd(M, R) = 1 is 
f*s selected. Division by R is preferably inexpensive, thus an optimal choice is R = 2 m if 

CO M ~ S"o m ' 2 ' * The m_residue of x is xR mod M - M ' = M_1 mod R is also computed. A 

[ft y 

tl function MRED(T) is provided that computes TR" mod M: This function computes the 

S ^ 15 normal representation of T, given that T is an m-residue. 

Cm 

r Method 1.3 MRED(T): computes a Montgomery reduction of T 
| T < RM, R = 2 m , M = Y!=l m i 2\ gcd(M , /?) = 1 

S 1. U = TM' mod R 

P 2. t = (T + UM) / R 

20 3. IF t>M RETURN t-M 

ELSE RETURN t 

The result of MRED(T) is t = TR" 1 mod M. 

Now to multiply two integers a and b in the transform domain, where their 
respective representations are (aR mod M) and (bR mod M), a product of the two 
25 representations is provided to MRED(T): 

MRED((aR mod M) • (bR mod M)) = abR 2 ^ 1 = abR mod M 
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For a modular exponentiation this step is repeated numerous times according to 
Method 1.1 or 1.2 to get the final result ZR mod M or P n R mod M. One of these values is 
provided to MRED(T) to get the result Z mod M or P n mod M. 

The initial transform step still requires costly modular reductions. To avoid the 
division involved, compute R 2 mod M using division. This step needs to be done only 
once for a given cryptosystem. To get a and b in the transform domain MRED(aR 2 mod 
M) and MRED(b*R 2 modM) are executed to get aR mod M and bR mod M. Obviously, 
any variable can be transformed in this manner. 

For a hardware implementation of Method 1.3: an m x m-bit multiplication and a 
2m-bit addition is used to compute step 2. The intermediate result can have as many as 
2m bits. Instead of computing U at once, one digit of an r-radix representation is 
computed at a time. Choosing a radix r, such that gcd(M, r) = 1 is preferred. Division by r 
is also preferably inexpensive, thus an optimal choice is r = 2 k . All variables are now 
represented in a basis-r representation. Another improvement is to include the 
multiplication A x B in the algorithm. 

Method 1.4 Montgomery Modular Multiplication for computing AB mod M, where 



M = J^&Ymt f ,m, e {o,1...2* - l}; B = 2^(2*)'*,,*, e {o,1...2* - l}; 



A,B < M;M <R = Z 1 ™ ; AT'= -M ~ l mod2* ; gcd(2* ,M) = 1 

1. S 0 = 0 

2. FOR i = 0 to m - 1 DO 

3. qi = (((Si + aiB) mod 2 k )M') mod 2 k 

4. S i+ i = (Si + qjM + ajB)/2 k 

5. END FOR 

6. IF S m > M RETURN S m - M 



The result of applying the method 1.4 is S m = ABR" 1 mod M. At most two k x k- 
bit multiplication operations and a k-bit addition is required to compute step 3 for a 



ELSE 



RETURN S, 
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radix 2 k . For step 4 two k xm- bit multiplication operations and two m + k— bit additions 
are needed. The maximal bit length of S is reduced to m+ k + 2 bits, compared to the 2m 
bits of Method 1.3. 

Method 1.5 is a simplification of Method 1.4 for radix r = 2. For the radix r = 2, 
the operations in step 3 of Method 1.4 are done modulo 2. The modulus M is odd due to 
the condition gcd(M, 2 k ) = 1. It follows immediately that M = 1 mod 2. Hence M' = -M" 1 
mod 2 also degenerates to M' = 1. Thus the multiplication by M* mod 2 in step 3 is 
optionally omitted. 

Method 1.5 Montgomery Modular Multiplication (Radix r = 2) for computing 
Montgomery Modular Multiplication for computing AB mod M, where 



M = I^G*)'",.*, e {0,1}; B = J£(2')%b, e {0,1}; A = j£(2')'a,.a l e {0,l}; 



A,B < M;M <R = 2 m ;gcd(2,M ) = 1 

1. S 0 = 0 

2. FOR i = 0 to m - 1 DO 

3. qj = (Sj + aiB) mod 2 

4. S|+i = (Sj + qjM + aiB)/2 

5. END FOR 

6. IF S m > M RETURN S m - M 



The final comparison and subtraction in step 6 of Method 1.5 would be costly to 
implement, as an m bit comparison is very slow and expensive in terms of resource 
usage. It would also make a pipelined execution of the algorithm impossible. It can easily 
be verified that S i+ i < 2M always holds if A, B < M. S m , however, can not be reused as 
input A or B for the next modular multiplication. If two more executions of the for loop 
are performed with a m+ i = 0 and inputs A, B < 2M, the inequality S m +2 < 2M is satisfied. 
Now, S m+ 2 can be used as input B for the next modular multiplication. 

To further reduce the complexity of Method 1.5, B is shifted up by one position, 
i.e., multiplied by two. This results in aj ■ B mod 2 = 0 and the addition in step 3 is 



ELSE 



RETURN S, 



m 
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avoided. In the update of Sj+i (Si + qM + a;B)/2 is replaced by (Sj + q f M)/2 + a;B. The 
cost of this simplification is one more execution of the loop with a m +2 = 0. The Method 
below comprises these optimizations. 

Method 1.6 Montgomery Modular Multiplication (Radix r=2) for computing A B mod 
5 M, where M = ££(2*) ' m n m i e B = XTo'* 2 ** e 

A = ^(2 k ) i a i ,a i e{0,l}; A,B<2M\M < /? = 2 ffl+2 ;gcd(2,M) =1 

1. S 0 = 0 

2. FOR i = 0 to m + 2 DO 

3. q; = (Si) mod 2 

10 4. S i+1 = (Si + qiM)/2 + ajB 

5. END FOR 



a 

CO The algorithm above calculates S m+3 = (2" (m+2) AB) mod M. To get the correct 

In 



■ 5- 



result an extra Montgomery modular multiplication by 2 2(m+2) mod M is performed. 
However, if further multiplication operations are required as in exponentiation 



U 1 2(m+2l 

%3 15 algorithms, it is better to pre-multiply all inputs by the factor 2 mod M. Thus every 
intermediate result carries a factor 2 m+2 . Montgomery multiplying the result by "1" 
eliminates this factor. 



The final Montgomery multiplication with "1" insures that a final result is smaller 
than M. 

20 High-radix Montgomery algorithm 

By avoiding costly comparison and subtraction operations of step 6 and changing 
the conditions to 4M < 2 km and A, B < 2M some optimisation results for implementing 
method 1.4 in hardware. The penalty is two more executions of the loop. The resulting 
method is as follows: 

25 Method 1.7 Montgomery Modular Multiplication for computing A*B mod M, 

whereM = J« (2 * ),w " w ." G f 0 ' 1 *" 2 * _1 ! 
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M = (M'mod2*)M,M = X^ 2 (2*)'m,,m, 6 f 0 ' 1 - 2 * 

< 2M;4M < 2* m ;M'= -M _1 mod 2* 
1. S 0 = 0 
5 2. FOR i = 0 to m - 1 DO 

3. qi = (Si + ajB) mod 2 k 

4. S,. +1 =(S,.+<7 ( .M+a,.fi)/2* 

5. END FOR 

The quotient qj determination complexity is further be reduced by replacing B by 
10 B • 2k. Since a;B mod 2 k = 0, step 3 is reduced to qj = Sj mod 2 k . The addition in step 3 is 
avoided at the cost of an additional iteration of the loop, to compensate for the extra 
CO factor 2 k in B. A Montgomery method optimized for hardware implementation is shown 

tV below: 

if*- 

u 

Cn Method 1.8 Montgomery Modular Multiplication for computing A • B mod M, where 
M = X m ; 0 3 (2*)'m,,m 1 . € {o,1...2* -l} 
M = (M'mod2*)M,M = J^U* >'^»*/ € f 0 ' 1 - 2 ' 

5 = ir 0 1 ( 2t) 'AAe {0.1...2* -l}; ^ = Xr =0 (2 t )'a ( .,« I € {o,1...2* -l},a m =0; 

A,B<2M;4M < 2*" ;M'=-M-' mod 2* 
1. So = 0 
20 2. FOR i = 0 to m - 1 DO 

3. qi = Si mod2 k 

4. S,. +1 =(S,.+<7 ( .M)/2*+a,.i? 

5. END FOR 

The final result is then Montgomery multiplied by 1 to eliminate the factors 
25 therein as discussed herein above. 



C3 
t0 



s 15 

□ 

in 
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In a thesis submitted to the Faculty of the Worcester Polytechnic Institute entitled 
Modular Exponentiation on Reconfigurable Hardware and submitted by Thomas Blum on 
April 8th, 1999 incorporated herein by reference, Thomas Blum proposed two different 
pipeline architectures for performing encryption functions using modular multiplication 
5 and Montgomery spaces: an area efficient architecture based on Method 1.6 and a speed 
efficient architecture. As target devices Xilinx XC4000 family devices were used. 

A general radix 2 systolic array uses m times m processing elements, where m is 
the number of bits of the modulus and each element processes a single bit. 2m modular 
multiplication operations can be processed simultaneously, featuring a throughput of one 
10 modular multiplication per clock cycle and a latency of 2m cycles. As this approach 
results in unrealistically large CLB counts for typical bit lengths required in modern 
m public-key schemes, only one row of processing elements was implemented. With this 
^3 approach two modular multiplication operations can be processed simultaneously and 
in the performance reduces to a throughput of two modular multiplication operations per 
[^15 2m cycles. The latency remains 2m cycles. 

\3 The second consideration was the choice of the radix r = 2 k . Increasing k reduces 

?3 the amount of steps to be executed in Method 1.8. Such an approach, however, requires 
j;0 more resources; The main expense lies in the computation of the 2 k multiples of M and B. 
■*y These are either pre-computed and stored in RAM or calculated by a multiplexer 
It 20 network. Clearly, the CLB count becomes smallest for r = 2, as no multiples of M or B 
have to be calculated or pre-computed. 

Using a radix r = 2, the equation according to Method 1.6 is computed. To further 
reduce the required number of CLBs the following measures are optionally taken: each 
unit processes more than a single bit. A single adder is used to precompute B+M and to 
25 perform the other addition operation during normal processing. Squares and 

multiplication operations are computed in parallel. This design is divided hierarchically 
into three levels. 

Processing Element Computes u bits of a modular multiplication. 
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Modular Multiplication An array of processing elements computes a modular 
multiplication. 

Modular Exponentiation Combine modular multiplication operations to a modular 
exponentiation according to Algorithm 1.2. 

5 Processing Elements 

Figure 10 shows the implementation of a processing element. 
In the processing elements the following registers are present: 

• M-Reg (u bits): storage of the modulus 

• B-Reg (u bits) : storage of the B multiplier 

O 

fg 10 • B+M-Reg (u bits): storage of the intermediate result B + M 

in 

• S-Reg (u + 1 bits): storage of the intermediate result (inclusive carry) 
%q • S-Reg-2 (u - 1 bits): storage of the intermediate result 

3 

i% • Control-Reg (3 bits): control of the multiplexers and clock enables 

kh 

C3 

v3 • ai,qj (2 bits): multiplier A, quotient Q 

15 • Result-Reg (u bits): storage of the result at the end of a multiplication 

The registers need a total of (6u + 5)/2 CLBs, the adder u/2 + 2 CLBs, the 
multiplexers 4 • u/2 CLBs, and the decoder 2 CLBs. The possibility of re-using registers 
for combinatorial logic allows some savings of CLBs. Mux B and Mux Res are implemented 
in the CLBs of B-Reg and Result-Reg, Muxi and Mux 2 partially in M-Reg and B+M- 
20 Reg. The resulting costs are approximately 3u + 4 CLBs per u-bit processing unit. That is 
3 to 4 CLBs per bit, depending on the unit size u. 

Before a unit can compute a modular multiplication, the system parameters have 
to be loaded. M is stored into M-Reg of the unit. At the beginning of a modular 
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multiplication, the operand B is loaded from either B-in or S-Reg, according to the select 
line of multiplexer B-Mux. The next step is to compute M + B once and store the result in 
the B+M-Reg. This operation needs two clock cycles, as the result is clocked into S-Reg 
first. The select lines of Muxi and Mux 2 are controlled by a t or the control word 
5 respectively. 

In the following 2(m + 2) cycles a modular multiplication is computed according 
to Method 1.6. Multiplexer Muxi selects one of its inputs 0, M, B, B + M to be fed in the 
adder according to the value of the binary variables a; and qj. Mux 2 feeds the u - 1 most 
significant bits of the previous result S-Reg 2 plus the least significant result bit of the next 
10 unit (division by two/shift right) into the second input of the adder. The result is stored in 
S-Reg for one cycle. The least significant bit goes into the unit to the right (division by 
two / shift right) and the carry to the unit to the left. In this cycle a second modular 
multiplication is calculated in the adder, with updated values of S-Reg 2 , &\ and q;. The 
second multiplication uses the same operand B but a different operand A. 



V3 the adder. This value is both stored into Result-Reg, as fed via S-Reg into B-Reg. The 
result of the second multiplication is fed into Result-Reg one cycle later. 

in 

P Figure 1 1 shows how the processing elements are connected to an array for 

O computing an m-bit modular multiplication. To perform the method for m bits with u bits 
r "20 processed per unit m/u + 1 units are used. Unito has only u - 1 B inputs as Bo is added to a 
shifted value Sj + qiM. The result bit S-Rego is always zero according to the properties of 
Montgomery's algorithm. Unitm/u processes the most significant bit of B and the 
temporary overflow of the intermediate result S i+ i. There is no M input into this unit. 



25 way. The control word, qi and aj are pumped from right to left through the units. The 

result is pumped from left to right. The carry-out signals are fed to the carry-in inputs to 
the right. Output SOOut is always connected to input S_0_ln of the unit to the right. 
This represents the division by 2 of the equation. 




At the end of a modular multiplication, S m+ 3 is valid for one cycle at the output of 



The inputs and outputs of the units are connected to each other in the following 
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At first the modulus M is fed into the units. To allow enough time for the signals 
to propagate to all the units, M is valid for two clock cycles. We use two M-Buses, the 
M-even-Bus connected to all even numbered units and the M-odd-Bus connected to all 
odd numbered units this approach allows to feed u bits to the units per clock cycle. Thus it 
takes m/u cycles to load the full modulus M. 

The operand B is loaded similarly. The signals are also valid for two clock cycles. 
After the operand B is loaded, the performance of the steps of Method 1.6 begins. 

Starting at the rightmost unit, unito, the control word, ai, and q\ are fed into their 
registers. The adder computes S-Reg-2 plus B, M, or B + M in one clock cycle according 
to a; and qj. The least significant bit of the result is read back as qi+i for the next 
computation. The resulting carry bit, the control word, a\ and qi are pumped into the unit 
to the left, where the same computation takes place in the next clock cycle. 

In such a systolic fashion the control word, ai, qi, and the carry bits are pumped 
from right to left through the whole unit array. The division by two in Method 1.6 leads 
also to a shift-right operation. The least significant bit of a unit's addition (So) is always 
fed back into the unit to the right. After a modular multiplication is completed, the results 
are pumped from left to right through the units and consecutively stored in RAM for 
further processing. 

A single processing element computes u bits of Sj+i = (Si + qj • M)/2 + aj • B. In 
clock cycle i, unito computes bits 0 ... u - 1 of Sj. In cycle i + 1, uniti uses the resulting 
carry and computes bits u . . . 2u - 1 of Sj. Unito uses the right shifted (division by 2) bit u 
of Si (S 0 ) to compute bits 0 ... u - 1 of S i+ i in clock cycle i + 2. Clock cycle i + 1 is 
unproductive in unito while waiting for the result of uniti. This inefficiency is avoided by 
computing squares and multiplication operations in parallel according to Method 1.2. 
Both p i+ i and z i+ i depend on z x . So, the intermediate result z x is stored in the B-Registers 
and fed with pi into the aj input of the units for squaring and multiplication. 

Figure 12 shows how the array of units is utilized for modular exponentiation. At 
the heart of the design is a finite state machine (FSM) with 17 states. An idle state, four 
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states for loading the system parameters, and four times three states for computing the 
modular exponentiation. The actual modular exponentiation is executed in four main 
states, pre-computationl, pre-computation2, computation, and post-computation. Each of 
these main states is subdivided in three sub-states, load-B, B+M, and calculate- 
5 multiplication. The control word fed into control-in is encoded according to the states. 
The FSM is clocked at half the clock rate. The same is true for loading and reading the 
RAM and DP RAM elements. This measure makes sure the maximal propagation time is 
in the units. Thus the minimal clock cycle time and the resulting speed of a modular 
exponentiation relates to the effective computation time in the units and not to the 
1 0 computation of overhead. 

Before a modular exponentiation is computed, the system parameters are loaded. 
The modulus M is read 2u bits at the time from I/O into M-Reg. Reading starts from low 
f p order bits to high order bits. M is fed from M-Reg u bits at the time alternatively to M- 
^ even-Bus and M-odd-Bus. The signals are valid two cycles at a time. The exponent E is 
M 15 read 16 bits at the time from I/O and stored into Exp-RAM. The first 16 bit wide word 

Cn 

tn 



o 

in 



from I/O specifies the length of the exponent in bits. Up to 64 following words contain 
the actual exponent. The pre-computation factor 2 2(m+2) mod M is read from I/O 2u bits at 
the time. It is stored into Prec-RAM. 



In state Pre-computel we read the X value from I/O, u bits per clock cycle, and 
M 20 store it into DP RAM Z. At the same time the pre-computation factor 22(m+2) mod M is 
read from Prec RAM and fed u bits per clock cycle alternatively via the B-even-Bus and 
B-odd-Bus to the B-registers of the units. In the next two clock cycles, B + M is 
calculated in the units. 

The initial values for Method 1.2 are available. Both values have to be multiplied 
25 by 2, which can be done in parallel as both multiplication operations use a common 

operand 2 2(m+2) mod M that is already stored in B. The time-division-multiplexing (TDM) 
unit reads X from DP RAM Z and multiplexes X and 1. After 2(m+3) clock cycles the 
low order bits of the result appear at Result-Out and are stored in DP RAM Z. The low 
order bits of the next result appear at Result-Out one cycle later and are stored in DP 
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RAM P. This process repeats for 2m cycles, until all digits of the two results are saved in 
DP RAM Z and DP RAM P. The result X • 2 m+2 mod M is also stored in the B-registers of 
the units. 

In state pre-compute2 the actual steps of Method 1.2 begin. For both calculations 
of Zl and PI ZO is used as an operand. This value is stored in the B-registers. The second 
operand ZO or PO respectively, is read from DP RAM Z and DP RAM P and "pumped" 
via TDM as a t into the units. After another 2(m + 3) clock cycles the low order bits of the 
result of Zl and PI appear at Result-Out. Zl is stored in DP RAM Z. PI is needed only if 
the first bit of the exponent eO is equal to "1". Depending on eO, PI is either stored in DP 
RAM P or discarded. 

In state compute the loop of method 1.2 is executed n - 1 times. Zj in DP RAM Z 
is updated after every cycle and "pumped" back as aj into the units. Pi in DP RAM P is 
updated only if the relevant bit of the exponent e\ is equal to "1". In this way always the 
last stored P is "pumped" back into the units. 

After the processing of e n .i, the FSM enters state post-compute. To eliminate the 
factor 2m-h2 from the result P n , a final Montgomery multiplication by 1 is computed. 
First the vector 0, 0, ... 0, 1 is fed alternatively via the B-even-Bus and B-odd-Bus into 
the B-registers of the units. P n is "pumped" from DP RAM P as ai into the units. After 
state post-compute is executed, u bits of the result P n = XE mod M are valid at the I/O 
port. Every two clock cycles another u bits appear at I/O. State pre-computel can be re- 
entered immediately now for the calculation of another X value. 

A full modular exponentiation is computed in 2(n + 2)(m + 4) clock cycles. That 
is the delay it takes from inserting the first u bits of X into the device until the first u 
result bits appear at the output. At that point, another X value can enter the device. With a 
additional latency of m/u clock cycles the last u bits appear on the output bus. 

Hereinbelow the function blocks in Figure 12 are explained. Figure 13 shows the 
design of DP RAM Z. An m/u x u bit DP RAM is at the heart of this unit. It has separate 
write (A) and read (DPRA) address inputs. The write-counter counting up to m/u 
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computes the write address (A). The write-counter starts counting (clock-enable) in sub- 
states B-load when the first u bits of Z x appear at data in. At the same time the enable 
signal of the DP RAM is active and data is stored in DP RAM. Terminal-count resets 
count-enable and write-enable of DP RAM when m/u is reached. The read-counter is 
enabled in the sub-states compute. When read-counter reaches its upper limit m + 2, 
terminal-count triggers the FSM to transit into sub-state B-load. The log2(m/u) most 
significant bits of the read-counter value (q out) address DPRA of the DP RAM. Every u 
cycles another value stored in the DP RAM is read. This value is loaded into the shift 
register when the log2(u) least significant bits of q out reach zero. The next u cycles u bits 
appear bit by bit at the serial output of the shift register. The last value of z\ is stored in a 
u-bit register. This measure allows us to select an m/uxu-bit DP RAM instead of an 
2m/uxu-bit DP RAM (m = 2x, x = 8, 9, 10). 

DP RAM P works almost the same way. It has an additional input ei, that activates 
the write-enable signal of the DP RAM in the case of ej = 1. 

Figure 14 shows the design of Exp RAM. In the first cycle of the load-exponent 
state, the first word is read from I/O and stored into the 10-bit register. Its value specifies 
the length of the exponent in bits. In the next cycles the exponent is read 16-bit at a time 
and stored in RAM. The storage address is computed by a 6-bit write counter. At the 
beginning of each compute state the 10-bit read counter is enabled. Its 6 most significant 
bits compute the memory address. Thus every 16th activation, a new value is read from 
RAM. This value is stored in the 16-bit shift-register at the same time when the 4 least 
significant bits of read counter are equal to zero. When read counter reaches the value 
specified in the 10-bit register, the terminate signal triggers the FSM to enter state post- 
compute. 

Figure 15 shows the design of Prec RAM. In state load-pre-f actor the pre- 
computation factor is read 2u bits at the time from I/O and stored in RAM. A counter that 
counts up to m/2u addresses the RAM. When all m/2u values are read, the terminal-count 
signal triggers the FSM to leave state load-pre-factor. 
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In state pre-computel the pre-computation factor is read from RAM and fed to 
the B-registers of the units. The counter is incremented each clock cycle and 2u bits are 
loaded in the 2u-bit register. From there u bits are fed on B-even-bus each positive edge 
of the clock. On the negative clock edge, u bits are fed on the B-odd-bus. 

A Speed Efficient Architecture 

The above design was optimized in terms of resource usage. Using a radix r = 2 k , 
k > 1, reduces the number of steps in Method 1.6 by a factor k. The computation of 
Method 1.8 is executed m + 3 times (i = 0 to m + 2) 

A speed efficient design is readily divided hierarchically into three levels. 

Processing Element Computes 4 bits of a modular multiplication. 

Modular Multiplication An array of processing elements computes a modular 
multiplication. 

Modular Exponentiation Combines modular multiplication operations to a modular 
exponentiation according to Method 12. 

Figure 16 shows the implementation of a processing element. 

The following elements are provided: 

• B-Reg (4 bits): storage of the B multiplier 

• B-Adder-Reg (5 bits): storage of multiples of B 

• S-Reg (4 bits): storage of the intermediate result Sj 

• Control-Reg (3 bits): control of the multiplexers and clock enables 

• ai-Reg (4 bits): multiplier A 

• qi-Reg (4 bits): quotient Q 
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Result-Reg (4 bits): storage of the result at the end of a multiplication 



B-Adder (4 bits): Adds B to the previously computed multiple of B 



B+M -Adder (4 bits) : Adds a multiple of M to a multiple of B 



S+B+M ~ -Adder (5 bits): Adds the intermediate result M ~ i 



S; tO B + 



B-RAM (16x4 bits): Stores 16 multiples of B 



• M ~ -RAM (16x4 bits): Stores 16 multiples of M 

The operation of the units is evident from the thesis of T. Blum, referenced above, 
and from a review of the diagrams. 

Figure 17 shows how the processing elements are connected to an array for 
computing a full size modular multiplication. 

Figure 18 shows how the array of units is utilized for modular exponentiation. 

Figure 19 shows the design of DP RAM Z. An mx4 bit DP RAM is at the heart of 
this unit. It has separate write (A) and read (DPR A) address inputs. Two counters that 
count up to m + 2 compute these addresses. The write-counter starts counting (clock- 
enable) in sub-states B-load when the first digit of Zi appears at data in. At the same time 
the enable signal of the DP RAM is active and data is stored in DP RAM. When m + 2 is 
reached, the terminal-count signal of the write-counter resets the two enable signals. The 
read-counter is enabled in sub-states compute. The data of DP RAM is addressed by q 
out of the read-counter and appears immediately at DPO. When read-counter reaches m + 
2, terminal-count triggers the FSM to transit into sub-state B-load. The last two values of 
zi are stored in a 4-bit register each. 

This measure allows us to choose a 100% utilized m x 4^-bit DP RAM instead of 
an only 50% utilized 2m x 4-bit DP RAM. DP RAM P works almost the same way. It 
has an additional input e i5 that activates the write-enable signal of the DP RAM in the 
case of ei = "1." 
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Since the above pipeline processor architectures embody many pipelined 
processing elements, it is often difficult and costly to synchronise each element to the 
clock source within a same integrated circuit. Therefore, the present invention is highly 
advantageous in reducing overall resource requirements by reducing clock distribution 
problems. Also, since in one direction addition is required while in the other direction 
multiplication is required, it is evident that more time is necessary along one path than the 
other and, so, time-averaging of the paths is possible in accordance with an embodiment 
of the invention. 

Numerous other embodiments may be envisaged without departing from the spirit 
or scope of the invention. 
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